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Abstract 

The calculations are performed for a broad range of the properties of Ca 
and Sr in the fee and bcc phases. A detailed information on the magnitude 
and character of temperature dependence of anharmonic effects in the lattice 
dynamics over the entire Brillouin zone (frequency shifts and phonon damping, 
Gruneisen parameters) is given. A detailed comparison of the computational 
results for the heat capacity and thermal expansion with the experimental 
data is carried out; the theoretical results are in good agreement with the 
experiment. 
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I. INTRODUCTION 



The investigation of anharmonic effects (AE) in lattice dynamics is a classical problem of 
solid-state physics. It is important, particularly, because of the role these effects can play in 
phenomena associated with structural phase transitions and melting in crystals (see, e.g.aH). 
At the same time, obtaining any information about the magnitude and scale of AEs from 
experiment and theory is a difficult problem. The experimental study of such "basic" AEs as 
the frequency shift and damping of phonons is very difficult and leads to a large uncertainty 
in the results (see, e.g., the data presented in@@ for bcc and fee metals, respectively). Up 
to now first-principles microscopic calculations of AEs have been performed for one point of 
the Brillouin zone (N) in the bcc phase of Zr and four points (N,P,ci;,G) in Moi. Detailed 
information about AEs in the entire Brillouin zone and their temperature dependence has 
been obtained in0i on the basis of pseudopotential theory for the bcc phases of alkali and 
alkaline-earth metals. For these metals the most striking manifestations of AEs are due to the 
"soft-mode behavior" (the anomalous temperature dependence of the phonon frequencies) of 
the S4 branch. It is of interest to calculate AEs for the "general position", i.e. for crystals 
which do not possess soft vibrational modes. Such crystals include most metals with close- 
packed structures, for example, fee. In a recent papeia the results of calculations performed 
for the anharmonic effects in lattice dynamics of Ir having the fee structure have been 
presented. In order to understand the specific features of AEs in lattice dynamics, related 
to structural phase transitions, it is interesting to investigate them in several phases for 
polymorphic metals. A classical example of structural phase transitions in metals is the 
fcc-bcc temperature transition in Ca and Sii. Therefore studies of AE features in these 
metals seem to be important. In the present paper the AE features in Ca and Sr lattice 
dynamics and thermodynamics have been investigated basing on microscopic calculations 
of the heat capacity, thermal expansion, temperature frequency shifts, damping of phonons 
and Gruneisen parameters. 



II. APPROXIMATIONS AND COMPUTATION PROCEDURE 

As noted in Introduction, at present, ab initio calculations of phonon frequency with 
regard to AE have been only done for some highly symmetric points of the Brillouin zone 
in a few metalsi. Consistent ab initio calculations of such quantities as phonon damping, 
as well as thermal expansion coefficients and other integral anharmonic effects in metals 
are now extremely difficult. Therefore in works reported so farBi'ioi particular models 
of interatomic interactions were used. For example, for alkali metalsi'^ and IiH the pair 
interaction approximation basing on pseudopotential theory was used, which gives a reliable 
description of a variety of lattice properties of these metals. Construction of a similar model 
for alkaline-earth metals appears to be difficult because of the specific features of electron 
structure. The fee phase of Ca and Sr has many van Hove singularities in the density of 
electron states near the Fermi level Ep and, therefore, for the adequate description of lattice 
properties allowance for related contributions to the total energy is important^. At the 
same time, the Fermi surface for these metals is close to that in the approximation of almost 
free electrons, and parameter \Vg\/EF (where g is the reciprocal lattice vector, \Vg\ is the 
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Fourier component of pseudopotential) is small (less than 0.1)0. As shown i n§, under these 
conditions the proximity of the Fermi surface to the Brillouin zone boundary essentially 
contributes to the elastic moduli and to frequencies of long-wave phonons with the wave 
vector 



|V^||-Ec — Ep 



q<qc = 9i\l (1) 

where Ec is the van- Hove singularity nearest to Ep, gi is the minimum vector of reciprocal 
lattice while the corresponding anomalous contributions to frequencies of phonons with q > 
Qc are much weaker. Therefore it may be expected that for the description of thermodynamics 
of fee phase of Ca and Sr for T > QnQc/gi ~ O.IGd (where 0_d is the Debye temperature) 
the ordinary approximation of pairwise interactions corresponding to the second order to 
perturbation theory for \Vg\/Ep would be apphcable. This model was successfully used in0 
for the description of the phonon spectra of Ca and Sr in the high temperature bcc phase. 

In the present paper the phonon spectra, thermodynamic and anharmonic properties 
were calculated using the model corresponding to the second order of perturbation theory 
for pseudopotential (see, e.g.SEl). For the latter the Animalou- Heine expression 



cos qvQ + U cos qvQ 

qro 



exp 



-0.03 



2k 



Fo 



(2) 



was used, where Z = 2 is the ion charge, e is the electron charge, rg, U are the pseu- 
dopotential parameters listed in table |, kp^ is the Fermi momentum at zero pressure. For 
the screening the Geldart-Taylor appmximat,o,Ji with the Ceperley-Alder expression for 
correlation energy was used (seeB for details). To eliminate the influence of fee band struc- 
ture peculiarities (the closeness of the Fermi surface to the bounadaries of the Brillouin 
zone), parameters Tq and U were fitted to the properties of high-temperature bcc phase: 
Q = flexpiT = Tg) and 6_d = 9^^(T = Tg) for Sr and Q = flexp(T = Tg) and C = C'^^^ for 
Ca. Here Tg is the temperature of bcc-fcc transition (see table |), index "exp" indicates the 
related experimental value, C = l/2(Cii — Cu) is one of the shear moduli. 

For the AE calculations the standard anharmonic perturbation theory was use dllii. With 
an accuracy up to terms of order, where x = (m/M)^/^ is the adiabatic parameter, m, 
M are the masses of electron and ion, respectively, the Hamiltonian of phonon subsystem 
can be written as: 

H = Ho + H,H + Hs + (3) 

where 

Ho = J2^>^Kh (4) 

A 

is the Hamiltonian of harmonic approximation, A = q^, q is the wave vector of phonon, ^ is 
the number of phonon branch, 6^, b\ are the phonon operators of creation and annihilation. 
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is the quasiharmonic Hamiltonian, Ui are the deformation parameters (in cubic crystals it is 
sufficient to allow for dilatation ui, dui = dlnQ, where Q is the lattice volume) 



H3= Yl V^'\Xl,X2,X3)QxrQx,Qxs 
A1A2A3 

Ai A2A3 



(6) 



H4 = ^ V^'-^'H^l, A2, A3, A4)Qai<5a2Qa3Qa4 
Ai A2A3A4 

= Yl ^^^\X,,X2,X3,X4)Ax,A,,A,^A,, 



Ai A2A3A4 

are the Hamiltonians of three- and four- phonon processes, respectively, 

1 



Q 



A 



Ax, Ax = bx + }P 



(7) 



(8) 



is the operator of phonon coordinate, —A = — q, ■C- Here and below we put the Planck 
constant as h = 1. Expressions for the amplitudes of three- and four-phonon processes 
${3)^ $(4)^ terms of the pseudopotential model used are shown in@. With an accuracy 
up to terms of order we have the following expressions for anharmonic shifts of phonon 
frequencies and damping (see Appendix): 



A, = Af ) + a1^) + a(^) 



A ^ ^A 



(9) 



A 



{qh) 



An 



-7a- 



n 



A 



(3) 



+ 



18 Y I^^'HA, -A2, -As)!' {(1 + ^A2 + N: 
A2A3 

+ {Nx, - Nx, 



V 



+ ^^A3 - ^\ 



V 



V 



^A3 - - ^^A - ^A. - ^^A 



(10) 



Af = 12 ^ <|.(^)(-A, A, A', -A')(l + 2iVv) 



(12) 



Vx = ISvr Y |*^'HA, -A2, -Aa)!' {(1 + Nx, + NxM^X2 + ^Xs - ^x) 
A2A3 

+ {Nx2 - Nx,) [6{iJX:i - WA2 - ^a) - 5(^A2 - ^A3 " ^a)]} 



(13) 



here 
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7A = (14) 

are the Gruneisen parameters, AQ is the change in the crystal volume due to the thermal 
expansion 

Nx = . , (15) 

is the Planck distribution function, V is the symbol of principle value. In the particular 
calculations the approximation 

V- ^ (16) 



6ix) ^ (17) 

was used for different small positive e with subsequent extrapolation e ^ 0. The details of 
the calculations, in particular, the summing up over the Brillouin zone is discussed in detail 
ini. 



III. CALCULATIONS OF PHONON SPECTRA IN HARMONIC 

APPROXIMATION 



The results of calculations of phonon spectra and Gruneisen parameters (|T^) for fee and 
bcc phases of Ca and Sr are shown in Figs. |l|-|^. As noted above, the theoretical model 
used is unapplicable for the description of phonon spectra in the fee phase of these metals 
in the nearest vicinity of point F because of unaccounted effects of Fermi surface proximity 
to the Brillouin zone boundaries0. It follows from Fig. ^ that qcC from Eq. (0) is of the 
order of 0.1 (71. Fig. ^ shows that the phonon spectrum in the bcc phase of Ca and Sr is 
qualitatively similar to phonon spectra in the bcc phase of alkali metalsi, exhibiting the same 
characteristic features: presence of "soft" S4 branch and prominent minimum for branch 
Fi. In the fee phase of these metals there are no soft mode anomalies. Fig. ^ displays such 
phonon spectra of bcc Sr for T = Tg with allowance for anharmonic frequency shifts and 
Corresponding experimental data in accordance withS. It is seen that there is a reasonable 
agreement between the calculated phonon spectra and experimental data. Unfortunately, 
any direct experimental data for uj{q) in the fee phase of Ca and Sr are unavailable. The 
comparison of theory with the experiment for this phase will be done in the next Section 
basing on the data on temperature dependence of the lattice heat capacity. As seen from 
Figs. H, ^ the Gruneisen model 7a =const is absolutely unapplicable, even qualitatively, for 
the description of bcc phases in alkali metals. In the fee phase the q-dependence is much 
weaker (provided that the vicinity of point F, where the model itself becomes unapplicable, 
is not considered). This seems to be because of lack of "soft modes" in the phonon spectra 
of fee phase. 
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IV. ANHARMONIC EFFECTS IN LATTICE DYNAMICS 



The calculation results for the anharmonic effects in the lattice dynamics for Ca and Sr 
are given in Figs. |6|-|T5|. The temperature dependencies of the frequency shift and phonon 



damping are shown for Ca as an example (Figs. ^ ^ [12], |T^); for Sr the anharmonic effects 
are somewhat stronger than for Ca, nevertheless all dependencies are similar. 
Figs. 1^ displays temperature dependencies of phonon frequency shifts: 

Acu(T) = uj{T)~uj{T = 0) 

= + Af ) + Af ^ (18) 

, (qh*) . (3*) . (4*) 



Here * , A)^ * , Aj^ * are values determined by Eqs. (|T0| - |13D respectively without the 



contribution of zero point oscillations, since such contributions enter both uj{T) and uj{T = 0) 
and are compensated in the difference (|T^). Note also that since fitting of the pseudopotential 
papameters was done at the volume Q = Qexp(T = Ts) (see Section Afl entering right 
hand side of Eq. ([TOD is determined by the following relation: AQ = ^1{T) — Q{Ts) 

It is seen that both in the bcc and in fee phases Au;(T) early (at T ^ 0.1 -r- 0.2Tpi) 
comes to high-temperature behavior Auj{T) ~ T. In the bcc phases of Ca and Sr the typical 
soft-mode behavior duj/dT > for S4 branch takes place (see the data for A^4 in Fig. |^). As 
pointed out inl3 the anharmonic effects in the bcc phases of alkaline-earth and alkaline metals 
are qualitatively similar in the whole, but, however, the scale of these effects in alkaline-earth 
metals is essentially larger. In the fee phase of Ca and Sr there are no "soft modes" and 
duj/dT < for all vibrational modes. 

Figs. 1^, ^ shows the q-dependence of the relative frequency shift in Ca and Sr. The 
comparison of the results for the fee phase (Fig. ^ with those obtained ini for Ir indicates 
that the positions of frequency shift minima and maxima in the Brillouin zone depends 
essentially on the ion charge; in It at Z = 4.5 maxima of S are reached at symmetric points 
X and L, while in alkaline-earth metals they are significantly shifted. Moreover, in contrast 
with the temperature dependencies S{T), there is a noticeable difference in the behavior of 
6{q) in Ca and Sr (cf., for example, the dependence in F — L direction in Fig. H). 

For the bcc phases of Ca and, particularly, Sr, 5(q) has a prominent maximum in the 
narrow region of q-space near point N for the soft branch S4 for which 6 > (see Fig. 

On example of Ca, Figs. |TU], |TT| show the contributions of three- and four- phonon 
processes to the frequency shift. In both phases these contributions have opposite signs in 
the larger part of the Brillouin zone (A3 < 0, A4 > 0) The mentioned above "soft mode" 
behavior of the branch E4 in the bcc phase is due to a sharp increase in A4 as in alkali 
metal^ and in BaTiO^-tjpe ferroelectricJil. 

In Figs. [T^, |TB| the temperature dependence of phonon damping in the symmetric points 
of Brillouin zone is shown for Ca as an example. First of all it should be pointed out that, 
in accordance with Eq. (|T3|) the damping does not vanish even at T = . It is due to 
three-phonon processes (phonon decay into two phonons) and has a smallness of order. 
It is seen from Figs. |T2|, that for Ca it is of the order of 0.2 -i- 0.5% of phonon frequency. 
As in the case of anharmonic frequency shift, it is seen that transition to asymptotic F ~ T 
is reached very early (at T ^ O.lTp;). 
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Figs. |14|, |T5| show the q-dependence of relative damping of phonons rj in Ca and Sr. In the 
fee phase rj has essentially nonmonotonous dependence on the wave vector, reaches values 
of about 25% in Sr and does not exceed 12% in Ca. In the bcc phase of Sr in the vicinity of 
point N f] reaches about 1 for the soft mode S4, which formally indicates inapplicability of 
the anharmonic perturbation theory in this narrow region. At the same time in the larger 
part of the Brillouin zone the anharmonic perturbation theory seems to be applicable. In 
view of this a problem arises on accurate separation of the contribution of the vicinity of 
point N to integral anharmonic characteristics such as free energy. This question will be 
considered elsewhere^. 



V. ANALYSIS OF EXPERIMENTAL DATA ON HEAT CAPACITY AND 

THERMAL EXPANSION 

A. Heat capacity 

The lattice heat capacity is determined by the phonon state density g{uj). Some char- 
acteristics of the latter, which can be used for comparing our description of phonon spectra 
with the experiment can be found from the analysis of temperature dependence of heat 
capacity at the constant pressure Cp{T). 

Here we carried out the separation of the lattice components from the heat capacity 
measured in the experiment and determined some average frequencies (momenta) of phonon 
spectrum, without using any model. 

In our analysis we used the following set of data: for Sr - direct experimental data 
from0, for Ca - the data from handbooks; in the region 2-50 K - the data from handbooklll 
which in the region 2-10 K is based onEl, and in the region 50-600 - the data taken from 
handbook^. 

In the analysis we neglected the temperature dependence of the electron heat capacity 
coefficient and assumed that the anharmonism is weak. The assumption on smallness of 
anharmonic contributions to the thermodynamical properties of metals is justified as it was 
checked by a direct calculation and found to be true even for alkali metals having soft phonon 
mo del. Under these assumptions the heat capacity measured experimentally at a constant 
pressure, Cp is described by the relations: 

Cp = Cph + Ca (19) 



Ca = {7 + {A- 7)(Cp,/3i?)2} J. (20) 

where 

Cph is the phonon component of the heat capacity in the harmonic approximation, 
Ca is the sum of contributions to the heat capacity, having the linear temperature de- 
pendence and coming from the anharmonic effects, thermal expansion of the lattice and 
conduction electrons; 

7 is the coefficient of electron heat capacity at low temperatures, 

A is the coefficient of linear in temperature term at high temperatures. 
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-R is the gas constant. 

The interpolation formula (^) for Ca gives the corresponding linear temperature asymp- 
totics both at low- and high-temperatures and ensures a smooth transition between the low 
and high temperature asymptotics in accordance with the law similar to Nernst-Linderman 
formula^. 

In the high-temperature region the Dhonon component of heat capacity was described 
by the expression proposed by NaumovE3 

C,. = SR{l--{^-^]\^('^y + J^]} (21) 




^ (1 - expiz)f V 12 240 y ^ ^ 

Here the asymptotic expansion of phonon heat capacity over small parameter z = Q/T is 
used: 



C,.^3flll-f;(!i^(^) ) (23) 



n=2 



where are the Bernulh numbers {B2 = 1/6, = -1/30, Bq6 = 1/42, Bg = -1/30, 
-Bio = 5/66 etc. For the odd n, beginning from n = 3, all -B„ = 0). 

In ([2T| ) the second and fourth order correction for Q/T are separately written, and func- 
tion v^(^) allows for all higher corrections in the " Einstein" approximation, i.e. assuming 
Qn = ^* for all n > 6. The values fi„ characterize momenta (mean frequencies) of phonon 
spectrum according to the relation: 

g{uj)uf'dujj J g{uj)duj (24) 

The value 7 entering Ca (pO|) was determined in the standard manner: by fitting of heat 
capacity in the low temperature region (particularly, 2-5 K for Ca and 5-15 K for Sr) by the 
relation Cp = 'jT + (3T^. The values 7, f3 as well as the limiting low-temperature value of 
the Debye temperature Blt related to f3 by relation [3 = 127r^i?/(5G|y) are listed in Table 
^ The estimations of 9^^ we obtained from the set of data used reasonably agree with 
those available in the literature (see01li). 

The value A together with ^2 , ^4 and fl^ were determined using the least square 



method by the fitting of heat capacity by relations (^, ^) in the temperature region 40-600 
K for Ca and 32-350 K for Sr. Within this temperature range these relations described 
the experimental results with least square deviation of the order of 1% for Ca and 0.3% for 
Sr. The parameters A, Q2, ^4 and determined by the least square method are shown 
in Table |1[ The table also gives the limiting high-temperature value of Debye temperature 
Qht = Qd(T) where Qlt >T <Tm (Tm is the melti ng t emperature) related to the second 
momentum of the phonon spectrum by: ^2 = Qht^/^J^^- 
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It should be pointed out that the lower value of the least square deviation for Sr is due 
to that over the whole temperature range the experimental data from the same source were 
used. On the other hand, its higher value for Ca seems to be due to non-smooth joining of 
the data at the transition from one source to another one. 

The analysis performed permitted the contribution of Cph in the harmonic approximation 
to be separated from the total heat capacity and the electron and anharmonic contributions 
to be left out. Some phonon spectrum momenta are expressed directly via phonon heat 
capacity integral; 



oo 



Cph 



<^>=2/ {^-^]dT (25) 



q fOO fl 



OO 

< uj-^ >= 0.138651 / ^ dT (27) 

OIL 



log cu>=^ % log ( T'^ dT (28) 
^ 7r2 Jo 3R ^ V 0-70702; ^ ^ 



We calculated these momenta with the integration over the experimental points in the 
temperature region 10-300 K and outside this region - over low temperature asymptotic 
Cph = l3T^ and high temperature asymptotic in the Debye spectrum model as inll. The 
mean frequencies corresponding to these momenta are listed in Table The meanings of 
fi_2) ^-1) ^1 corresponding to that of Qn in relation ^ for n = -2, - 1 and 1 while the 
meaning of Qiog is determined by relation: 

\og{niog) = — — — r- — 

oo oo 

,(.)log.;^^ / Jgju^^^ ^29) 



00 UJ 



The contribution of Cph in the harmonic approximation separated from the total heat ca- 
pacity is shown in Figs. |T6|, IT7. 



B. Thermal expansion 

The coefficient of thermal expansion a{T) contains two contributions: the electron one 
Ue and the lattice one ai. In the low temperature region the electron contribution is linearly 
dependent on temperature, while the lattice contribution is proportional to T^. This makes 
it possible to separate these contributions fitting the experimental data by the dependence: 
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a(T) = DT + ET^ (30) 

In the papeiii the electron contribution to thermal expansion of Ca and Sr is determined. 
For the determination of the lattice contribution to thermal expansion we used the following 
set of experimental data: in the low temperature region — the data from0 and in the 
high temperature region — the data fromi3 for Ca and@ for Sr. The phonon contribution to 
thermal expansion was determined as the difference of the experimental data and the electron 
contribution. It was assumed that the electron contribution is linear by temperature over 
the whole temperature region, its value being determined by low temperature asymptotic 
in@. The lattice contribution to thermal expansion separated from the experimental data 



is shown in Figs. IS, IS. 



The lattice contributions to the heat capacity and the coefficient of thermal expansion, 
separated from the experimental data were used for the determination of the Gruneisen 
lattice parameter by the formula: 

G = (31) 

where B is the bulk modulus, a is the temperature coefficient of linear expansion, Vm is the 
molar volume. The temperature dependence of the Gruneisen parameter 7(T) is shown in 
Figs. ITR ra. 



VI. CALCULATIONS OF HEAT CAPACITY AND THERMAL EXPANSION 



Figs. |1^, |1^ give the computation results of lattice heat capacity of Ca and Sr in the 
harmonic approximation and the effective Debye temperature. It is seen that, with the 
exception of the region of lowest temperatures where the phonon spectrum model itself be- 
comes inapplicable (see Section |I|), the theory is in a perfect agreement with the experiment. 
It is especially well seen in Table ||, where the computation data are compared with the 
momenta of phonon density of states, found from the experimental data. 

Figs. |18|, [1^ displays the computation results for the temperature dependence of the 
volume Ai7(T) and of the Gruneisen parameter 7(T) for Ca and Sr. It is related to the 
Gruneisen parameters by the following expression: 

7(T) = (32) 

5^iV,(l + iV,)^2 

A 

One can see that the temperature dependence of the macroscopic Gruneisen parameter is a 
rather smooth. At the same time, our calculations show a strong temperature dependence 
of microscopic Gruneisen parameters for the soft branch in BCC phases of Ca and Sr (see 
Table |IT|). It would be very interesting to check this prediction experimentally. 
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VII. CONCLUSIONS 



In the present paper a simple, and at the same time, adequate model of the interatomic 
interactions is proposed, which describes with a reasonably high accuracy the lattice proper- 
ties of fee and bcc phases of Ca and Sr. As a result, the dependence of AE in the dynamics 
of lattice on its geometry could be investigated. The comparison of the results obtained for 
the fee phase of Ca and Sr with earlier results for Iri demonstrates an essential dependence 
of AE on the character of itneratomic interactions for the given crystal structure. As the 
question on the AE scale, particularly, for metals with polymorphic transformation (in this 
case bcc-fcc) is one of the key problems in the lattice dynamics, performance of further 
experimental investigations is desirable. 



APPENDIX: 

Here we shall give a short derivation of equations for the temperature shift of frequencies 
and phonon damping (PHT3|). In contrast to the basing wor and subsequent publications, 
where the diagram method for Matsubara Green functions was used, we employ the method 
of two time Green function£ll. We shall begin with the chain of equations of motion for 
delayed two time Green functions 

< A|5 < > (Al) 

where 

oo 

<A|fi>^=-i j dtexp[i{uj + ie)t]{[A{t),B{0)]) \^^^^ (A2) 



The bracketed expression indicates the commutator of operators, and that within the French 
quats — the averaging over the Gibbs ensemble. For the Hamiltonian we have the 

accurate equation: 

{uj - ujx) < bx\b1^ >^ =1 + 3J2 A2, A3) < Ax.Ax.lb-^ >^ 

A2A3 

+4 J2 A2, A3, A4) « Ax,Ax,Ax,\bt »^ (A3) 

A2A3A4 

With an accuracy up to the Green function describing the contribution of four-phonon 
processes to (|A3|) can be decoupled 



< Ax,Ax,Ax,\b+ >^ = {Ax,Ax,) < Ax^lb^ >^ +{Ax,Ax,) < Ax,\b+ >^ 

+ {Ax,Ax,) <^Ax,\b+ >^ 
= SxiX-si^ + 2Nx2)5x4X + '^A2A4(1 + 2NxJ5x3X 

+6x,x,{l + 2Nx,)6x,x+ < bx\b^ >^ (A4) 
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This is immediately followed by the result (|T^) for the contribution of four-phonon processes 
to the frequency shift. For the determination of the contribution of three-phonon processes 
we must write the chain of equations for the Green-functions entering the following relation 

+ < >^ + < b-x^b-xM >^ (A5) 

and make there decoupling. Upon separation of the real and imaginary parts of Green 
function -C &a|&a obtain immediately the result (|9|-[T3|). The quasiharmonic frequency 

shift proves to be then as expressed through the amplitudes of three-phonon processes: 

Af ) = -36 V i±^$(A, -A, X^mXs, -A3, -A^) (A6) 

As far as we know this expression is new, although the relation itself of thermal expansion 
to three-phonon anharmonism was qualitatively discussed as early as in the classical book 
by Paierl£i. 
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FIGURES 



FIG. 1. Phonon spectrum of fee phase of Ca (solid line) and Sr (dashed line) in the harmonic 
approximation in units of ion plasma frequency ojpi = {AirZe^ /Miiy^'^ . Indices 1 and 2 denote 
longitudinal and transverse branches respectively. 

FIG. 2. Phonon spectra for bcc phase of Ca (solid line) and Sr (dashed line) in units ujpi in the 
harmonic approximation. 

FIG. 3. Phonon spectra for bcc phase of Sr at T = T^; solid line is the calculation (with the 
account for anharmonic shift, see below), points - the experimental data fromlll. 



FIG. 4. Gruneisen parameters for fee phase of Ca (solid line) and Sr (dashed line). 

FIG. 5. Gruneisen parameters for bcc phase of Ca (solid line) and Sr (dashed line). 

FIG. 6. Temperature dependence of frequency shift determined by Eq. (^) at symmetric 
points of the Brillouin zone for fee Ca; copi = {ATiZe^ /Mfl)^/'^ is the plasma ion frequency, Tpi is 
the corresponding temperature (for Ca Tpi=476 K.) The symbols correspond to different vibration 
modes: o ~ Wi, A — W2, D — Xi, — X2, — Li, ★ — Lg. 

FIG. 7. Temperature dependence of frequency shift at symmetric points of the Brillouin zone 
for bcc Ca; The symbols correspond to different vibration modes: o — Ni, (} — N3 A — A'4, □ 
— H,'^ — P. 

FIG. 8. Change of the relative frequency shift 6^i^{T) = Auj^i^{T)/uj^i^{0) in the Brillouin zone 
for fee phases of Ca and Sr at the temperature of structure transition T = Tg (Tg = 726 K for Ca 
and 930 K for Sr), Atj^k(T)/ is determined by Eq. (0). Solid line — Ca, dashed line — Sr. 

FIG. 9. The same as Fig. |8| for bcc phase of Ca and Sr. 
FIG. 10. Change in the Brillouin zone of the anharmonic frequency shift determined by formu- 



lae dnl), dll) for fee Ca at T = T^. 

FIG. 11. The same as Fig. |l^ for bcc phase of Ca. 

FIG. 12. Temperature dependence of phonon damping at symmetric points of Brillouin zone 
for fee Ca. Tpi and the symbols are the same as in Fig. 0. 
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FIG. 13. The same as Fig. for bcc phase of Ca. The symbols are the same as in Fig. |^ 



FIG. 14. Change in the Brillouin zone of relative damping r]^\i{T) = r^k(^)/'-^5k(0) for fee Ca 
and Sr at T = Tg. Solid line — Ca, dashed line — Sr. 



FIG. 15. The same as Fig. 14 for bcc phase of Ca and Sr. 



FIG. 16. Phonon contribution to the lattice heat capacity of Ca; in the insert the temperature 
dependence of effective Debye temperature is shown. is the result of experimental treatment (see 
the text); solid line — the calculation result for the bcc phase, dashed line — for the fee phase. 



FIG. 17. The same as Fig. 16 for Sr 



FIG. 18. Phonon contribution to the thermal expansion of Ca; in the insert the temperature 
dependence of Gruneisen parameter is shown. The meaning of symbols and lines is the same as in 
Fig. H. 



FIG. 19. The same as Fig. 18 for Sr. 
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TABLES 



TABLE I. Model parameters. 



Metal 


Ts 




copi (THz) 


Tm (K) 


ro (au) 


U 


Z 


Ca 


726 


303.4 


9.91 


1112 


2.545 


0.636 


2 


Sr 


930 


394.6 


5.88 


1045 


3.024 


0.787 


2 



TABLE II. Phonon spectra momenta, 
momenta Ca Sr 





experiment 


theory 


theory 


experiment 


theory 


theory 




FCC 


FCC 


BCC 


FCC 


FCC 


BCC 


^log, K 


136 


138 


128 


81 


78 


67 




130 


142 


121 


77 


77 


64 


K 


152 


158 


148 


89 


88 


80 


ni, K 


181 


185 


182 


102 


102 


99 


^2, K 


195 


193 


193 


107 


107 


106 


^4, K 


214 


207 


208 


116 


114 


114 


n„ K 


220 






120 






mJ 
'' molK^ 
a mJ 

Olt, K 


3 

1.87x10-2 
232 






13.8 

8.9 xlO-2 
138 






A mJ 

molK'^ 
@HT,K 


9.1 

252 


249 


249 


6.9 

139 


138 


137 



TABLE III. Temperature dependence of Gruneisen parameter. 



A^i N3 N4 R P 

Ca, BCC T=726 K 

70 1.60 1.12 3.60 1.31 1.66 

A7(T) 0.059 0.12 -2.87 0.24 -0.13 

7(T) 1.66 1.24 0.73 1.55 1.53 

Sr, BCC T=930 K 

70 1.69 1.16 5.44 1.37 1.77 

A7(r) 0.13 0.31 -20.0 0.33 -0.37 

7(r) 1.82 1.47 -14.6 1.7 1.40 
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